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The relativistic nature of the electron motion underlies the intrinsic part of the anomalous Hall 
effect, believed to dominate in ferromagnetic (Ga,Mn)As. In this paper, we concentrate on the 
crystal band structure as an important facet to the description of this phenomenon. Using different 
k.p and tight-binding computational schemes, we capture the strong effect of the bulk inversion 
asymmetry on the Berry curvature and the anomalous Hall conductivity. At the same time, we find 
it not to affect other important characteristics of (Ga,Mn)As, namely the Curie temperature and 
uniaxial anisotropy fields. Our results extend the established theories of the anomalous Hall effect 
in ferromagnetic semiconductors and shed new light on its puzzling nature. 



I. INTRODUCTION 

The anomalous Hall effect was first observed in fer- 
romagnets by Hall himself.— Next to the usual Lorentz 
term, a voltage proportional to magnetization appeared 
- much smaller, but still too large to be explained by 
an internal magnetic field. Over the following years, this 
ostensibly plain dependence was to unveil the whole cor- 
nucopia of physical phenomena, all arising from the rel- 
ativistic coupling of the charge and spin current A 
class of them is related to the spin asymmetry of car- 
rier scattering, viz. the skew-scattering^— and the side 
jump processor— Also higher-order effects in the scat- 
tering amplitude were predicted.— & 10 ' 1 1 However, it is 
the "intrinsic" mechanism, first proposed by Karplus and 
Luttinger^ that is believed to play a key role in the Hall 
effect of ferromagnetic semiconductors J^r— Quite unusu- 
ally, it does not concern the changing of the occupations 
of Bloch bands by scattering on impurities. Rather, it re- 
sults from the interband coherence caused by the univer- 
sal tendency of physical systems to progressively increase 
the indeterminacy of their state. Hence, the topological 
theory of the Berry phased was found to provide an am- 
ple description ! 20 ' 21 

The topological model of the intrinsic anomalous Hall 
effect describes the linear response of the carrier Bloch 
function ip to the applied electric field E. It consists in 
the drift of the k vector in the reciprocal space, during 
which tpCk) acquires a geometrical phase factor in addi- 
tion to the dynamical one. The Berry phase in the first 
factor can be expressed as the action of the vector po- 
tential, with its curl called the Berry curvature, fi(k). 
The latter is a well-defined gauge-invariant quantity of- 
ten pictured as a non-homogeneous magnetic field living 
in fc-space. It produces a local equivalent of the Lorentz 
force, the so-called anomalous velocity term, — eE x fl, in 
the semiclassical equations of motioni 12,20 This term con- 



tributes to the stationary part of the Boltzmann trans- 
port equation (hence, it does not depend on the trans- 
port relaxation time), producing a dissipationless current 
transverse to E. The anomalous Hall conductivity of this 
current is proportional to the ensemble average of the car- 
rier Berry curvature, (f2). Since ft depicts the changes of 
the spin polarization during the carrier transport by the 
electric field, which are caused by the spin-orbit coupling, 
it changes sign under time-reversal symmetry. Thus, to 
obtain a finite value of (f2), this symmetry of the system 
must be broken. 

The above semiclassical approach, taking into account 
the complete geometrical Bloch state description, leads 
to an intuitive picture of the origin and mechanism of the 
intrinsic anomalous Hall effect in ferromagnetic semicon- 
ductors. It applies to the weak scattering regime, where 
it was proven to be formally equivalent to more sys- 
tematic quantum-mechanical techniques^ In this frame- 
work, the anomalous Hall conductivity was calculated for 
the p-d Zener model^ 3 . by Jungwirth et alJ^ The band 
structure of this model is composed of the six hole bands 
described by the Kohn-Luttinger Hamiltonian with the 
mean-field spin splitting, neglecting the spin-orbit in- 
duced Rashba (linear in fc) and Dressclhaus (fc 3 ) terms. 
While the former, together with all terms linear in fc, 
does not generate the spin current^! the latter does^ 
which has not been hitherto studied in (Ga,Mn)As and 
related ferromagnets. 

In this paper, we investigate the intrinsic anomalous 
Hall effect (for our purposes called the AHE) in the p-d 
Zener model of a diluted ferromagnetic semiconductor, 
focusing on (Ga,Mn)As. This problem requires a com- 
plete description of the host band structure. We demon- 
strate it numerically by using different k.p and tight- 
binding computational schemes described in Sec. [Til We 
employ the 6-ban d 23 ' 26 ~ — and the 8-band k.p model in- 
cluding the Dresselhaus splitting ) 29 ' 30 and two empirical 
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tight-binding paramctcrization s 31 ' 32 (spds* and sps*) to 
describe the GaAs band structure. The Mn Ga substitu- 
tions are introduced within the mean-field and virtual- 
crystal approximations. In Section IlIIl we first calibrate 
the models so as to obtain the agreement of two impor- 
tant characteristics of (Ga,Mn)As, the Curie tempera- 
ture and uniaxial anisotropy. Then, we calculate the 
Berry curvature and the anomalous Hall conductivity to 
reveal qualitative differences between the calibrated mod- 
els. We even report a negative conductivity sign within 
the new approaches, which was not observed in the pre- 
viously employed 6-band k.p models This result relates 
to the inversion asymmetry of the zinc-blende lattice, in- 
herited by the Berry curvature. We provide the physical 
interpretation of our findings and briefly discuss their ex- 
perimental implications. 

II. THEORETICAL APPROACH 

We investigate AHE in a Hall sample of ferromagnetic 
(Ga,Mn)As with the electric field E || x and the magnetic 
field applied along the [001] direction. This setup yields 
the anomalous conductivity 

e 2 

<*xy = -y^{^z) , (1) 

where (tt z ) = J2k,n ^z( n > k )/™,k, and /„, k is the Fermi- 
Dirac distribution associated with the band n and wave 
vector k. The positive values of a xy mean that the 
anomalous Hall voltage has the same sign as in the ordi- 
nary Hall effect. 

The Berry curvature in (p} is given by 

fi z (n,k) = 2lui{d ky u n ^\d kw u n ^) , (2) 

or by the equivalent Kubo formula (derived by differenti- 
ating the Schrodinger equation over k, which makes sense 
in our finite-dimensional Hilbcrt space) 

where c nn i = (it n ,k|Sk-ffk|un',k), and w raj ic are the peri- 
odic parts of the Bloch states with energies e n ^. For- 
mula ((21) may carry large error even when we describe 
the Bloch wave functions quite accurately, because it in- 
volves their derivatives. For instance, in the 6-band k.p 
model , 27 ' 28 we obtain an almost perfect description of the 
p-typc bands around the T point, but their derivatives in 
general include significant contributions from the states 
outside this space. On the other hand, the sum in ([3]) 
goes over all bands, not just the hole p-type ones. Even 
the detailed description of these bands only is, therefore, 
not sufficient to calculate the Berry curvature accurately. 
The model used must also have enough room to allow for 
the inversion symmetry breaking, an important property 
of GaAs latticed 



For the above reasons, we expect the multiband tight- 
binding model s 31 ' 32 of the host band structure to be the 
most appropriate for the description of the Berry cur- 
vature. They automatically take into account the lack 
of inversion symmetry, as they distinguish Ga and As 
atoms. Contrary to the perturbativc k.p methods, they 
properly describe the Bloch states away from the center 
of the Brillouin zone, which makes them better suited to 
high hole concentrations. We use the spds* Jancu 3 ^ and 
sps* di Carlo^ 2 - parameterizations, basing our numerical 
tight-binding implementation on the code by Strahberger 
et al.^ employed previously in the studies of spin trans- 
port properties in modulated (Ga,Mn)As structures! 35 ' 36 

The impact of the inversion symmetry breaking on 
the Berry curvature is additionally tested in the 8-band 
k.p model with the Dresselhaus term included, following 
Ostromeki^ He found that the magnitude of the Dressel- 
haus spin splitting of the conduction band in GaAs de- 
pends on the values of two poorly known parameters A' 
and B describing the spin-independent and spin-orbit re- 
lated k.p interaction of the conduction band with remote 
bands, respectively. Two sets of A' and B values repro- 
duced the experimental magnitudes of spin splittings^ 
We have adopted the set for which A' = instead of 
the alternative one with A' = — 14.7eVA, which ap- 
pears unrealistic^ 7 - For the chosen parametrization, at 
S = (inversion symmetry preserved) the remaining pa- 
rameters correspond to the standard 6-band k.p model: 
7i = 6.85, 72 = 2.9, 73 = 2.1, and spin-orbit splitting 
A so = 0.341 eV. The fc-dependent part of the spin-orbit 
interaction has negligible effect on the investigated quan- 
tities, so we neglect it for clarity. 

The dispersion relations of the top of the valence band 
calculated by the above methods are compared in Fig.Q] 
There is a very good agreement between the most popular 
6-band k.p and the most detailed spds* tight-binding, as 
well as the 8-band k.p model, while the results obtained 
within sps* parameterization differ slightly. 

The biaxial strain is included in particular computa- 
tional schemes by adding an appropriate Bir-Pikus ma- 
trix to the 6-ban d 27 ' 28 and to the 8-band k.p model* 3 - 8 - 
or by changing the atoms arrangement in the tight- 
binding approaches, according to the strain tensor val- 
ues: e xx = E yy = Aa/a and e zz = —Icyxj c\\s XXy where 
Aa is the strain-induced change of the lattice constant 
a, and cnjcw = 0.453 is the ratio of clastic moduli. 
Additionally, the on-site energies of the d orbitals in the 
spds* parameterization depend linearly on the strain ten- 
sor values i^L 

The AHE current flows in ferromagnetic (Ga,Mn)As, 
in which a part of Ga atoms is substituted by Mn ions 
producing strong p-d hybridization. Each of them si- 
multaneously forms a magnetic moment 5 = 5/2 and 
creates one moderately bound hole. The holes fill the 
GaAs valence band, which is mostly of As p character, 
from top to bottom. Their concentration p is taken 
into account by adjusting the Fermi energy Ep. To 
do it efficiently, we assume that the crystal has a fi- 
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FIG. 1: [color online] The top of the GaAs valence band with 
spin splitting A = —0.1 eV calculated by different methods. 
Horizontal lines denote the positions of the Fermi level Ef for 
the hole densities 0.3 and 1.2 nm -3 in the spds* tight-binding 
model. 



III. RESULTS 
A. Curie temperature and uniaxial anisotropy 

We begin with the comparison of the band structure 
models, looking at how they describe the two important 
characteristics of the p-type hole bands, the Curie tem- 
perature Tc and uniaxial anisotropy field H nn . 

The four models employed provide similar values of 
the Curie temperature Tc, presented in Fig. [2] for three 
different Mn contents function of the hole concen- 

tration p. While the 6-band k.p and spds* tight-binding 
model give virtually identical Tc values, the remaining 
ones exhibit some differences due to their parameteriza- 
tion flaws (Fig.[T|). The slight discrepancy between the 
two best results for high hole concentrations is resolved 
in favour of the more universal tight-binding approach. 
Since Tc is proportional to the thermodynamic spin den- 
sity of states ; 23 ' 27 we conclude that a mutually consis- 
tent description of the relevant valence bands is achieved 
throughout. 



nitc and very small volume V ~ 10 7 ciq, and find the 
minimum of the first n = pV occupied hole states' en- 
ergies. The holes mediate ferromagnetic order between 
Mn spins via exchange interactions. Within the mean- 
field and virtual crystal approximations, which we em- 
ploy in all computational schemes, these interactions are 
spatially averaged. The constant molecular field of Mn 
spins creates a fc-dependent Zceman-likc splitting of the 
host bands (Fig.QJ. This splitting for heavy holes in 
the T point is given by A = xN n (S)/3, where x is the 
part of cation sites iV substituted by Mn with an av- 
erage spin (S), and /3 = —0.054 eVnm 3 is the p-d ex- 
change integral^ Usually, the A value is smaller than it 
would result from the nominal Mn concentration. With 
increasing Mn doping, a part of Mn atoms occupies in- 
terstitial positions and tends to form pairs with substi- 
tutional ones, characterized by very small net magnetic 
moment^ They can be removed by annealing, but the 
effective Mn concentration will remain smaller. Addition- 
ally, unintentional defects such as Mn interstitials and As 
antisites are double donors and reduce the hole concen- 
tration. 

In both k.p and tight-binding models, we have been 
able to compute the derivatives of the Hamiltonian ma- 
trix in ([3]) analytically, which significantly improves the 
accuracy of our results. Formula ([2]) is equally suitable 
for numerical computation, if we overcome the problems 
created by the wavefunction phase gauge freedom, which 
is cancelled analytically, but not numerically. One should 
simply fix the phases of the relevant wave functions by 
dividing each of them by the phase factor of its first non- 
zero basis coefficient. 
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FIG. 2: [color online] Computed magnitudes of Curie tem- 
peratures as a function of the hole concentration according 
to various band structure models for Gai-^Mn^ As with the 
Mn content x = 2%, 5% and 8%. Vertical lines indicate the 
maximal experimentally realisable hole densities for given x. 
On this and the following graphs, the numerical data were 
generated more densely than the markers, which are visual 
aids only. 

A non-trivial comparison between the used models is 
provided by evaluating the magnitude of the uniaxial 
magnetic anisotropy field H un brought about by biax- 
ial strain i 27 ' 40 This anisotropy is driven by the presence 
of the spin-orbit interaction in the carrier band. We cal- 
culate the magnitude of H un as proportional to the dif- 
ference of the total carrier energy for the easy and hard 
magnetization directions under 1% tensile or compressive 
strain (see Ref.l27h. As presented in Fig.[3l in a region of 
intermediate hole concentrations, the easy axis takes the 
[001] direction for tensile strain, while it is in the (001) 
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plane for compressive strain. The situation is opposite 
for lower and higher hole concentrations. These results 
agree between the models, especially for the spds* tight- 
binding and the 6-band k.p calculations. Consequently, 
all models handle similarly well the spin-orbit splitting 
of the p-type valence bands. 
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FIG. 3: [color online] The amplitude of the uniaxial 
anisotropy field H^n (divided by saturation magnetization, see 
Ref.l27h for compressive (e xx — —1%) and tensile (e xx = 1%) 
biaxial strain in (Ga,Mn)As with spin splitting A = — 0.15eV. 
The effect of shape anisotropy is neglected. In the central 
range of hole concentrations (II), the easy axis is in-plane and 
perpendicular to the plane for compressive and tensile strain, 
respectively, while the reorientation transition is expected in 
either low or high concentration regimes (I or III). 



B. Berry curvature and the related conductivity 

According to the previous chapter we can character- 
ize both Tc and H un as static quantities, which depend 
on the properties of the six occupied p-type bands only. 
Details of the other bands' structure, in particular the 
Dresselhaus k 3 splitting, do not influence their values. In 
marked contrast, the derivatives and interband elements 
in the Berry curvature formulas, (0 and ([3]), express the 
dynamic character (related to the carrier drift caused by 
electric field) of the AHE and lead to qualitative differ- 
ences between the models. Below we demonstrate the 
effect of the bulk inversion asymmetry on the Berry cur- 
vature and consequently on the anomalous conductivity 
trends. 

The 6-band k.p model describes the diamond lattice 
structure. Since the Kohn-Luttinger Hamiltonian it uses 
is invariant under space inversion, which is unitary, we 
have f2(— k) = fi(k). On the other hand, the antiunitar- 
ity of time reversal operator leads to fi(— k) = — S7(k) in 
the presence of the corresponding symmetry. Thus, the 
Berry curvature in this model is always symmetric and 
vanishes in the absence of magnetic fields, as presented 
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FIG. 4: [color online] Berry curvature of majority heavy and 
light hole bands in (Ga,Mn)As calculated using a) the 6-band 
k.p, b) 8-band k.p and c) spds* tight-binding model, with and 
without the spin splitting A. Inset: Divergences in the Berry 
curvature. 



The 8-band k.p model contains, in addition, the s-typc 
conduction band with the Dresselhaus spin splitting in- 
cluded by the use of Lowdin perturbation calculus . 29 ' 30 
It results from the inversion symmetry breaking in the 
zinc-blende structure^ 3 , and thus, in the presence of the 
time-reversal symmetry, leads to non-vanishing antisym- 
metric Berry curvatures (Fig.0t>). Then, the related k 3 
energy term in the conduction band spectrum accounts 
for a non-zero spin current within the intrinsic spin Hall 
effect 3^ When the magnetic field is on, the significant 
asymmetry of the curvatures can still be observed. 

The multiband tight-binding methods give us the de- 
tailed band parameterization and introduce realistic sym- 
metries of the crystal lattice in a natural way. Figure 2J; 
presents the Berry curvatures obtained using the spds* 
parameterization. Their symmetry is similar to the 8- 
band k.p model, but the shape differs (especially for 
A = 0), pointing to the sensitivity of the spin topological 
effects to the subtleties of the band structure. 

An interesting effect is the formation of so-called dia- 
bolic points corresponding to the energy bands' crossings, 
best visible for fi a (k) in the k || [001] direction (Fig.gt, 
inset). A commonly held view is that it is them which 
are the source of the anomalous Hall conductivity. Even 
though the degeneracies of states in the fc-space do pro- 
duce a nontrivial Berry potential, it is easy to show that 
their contributions to a xy vanish for T — 5- K. The coef- 
ficients c x ^ n , in ([3]) are the matrix elements of Hcrmitian 
operators dk x Hk> hence d^ n , = (c£{ n )*. The conductiv- 
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ity is thus proportional to the sum 

2 im E E (f cL f:' n ^ - • (4) 

k n<n' ' ' 

For e Mj k 7^ £n'.k, a component of the above sum with 
given (k, n, n') has a non-zero contribution to a xy only if 
fn.k =/= /n',kj which for T — ► happens when one state is 
above and another below the Fermi level Ep. (Exploit- 
ing this observation in numerical computations ensures 
fast convergence of calculated integrals.) The component 
corresponding to the bands' crossing is thus zero, since 
fnM = fn',k m a neighbourhood of the diabolic point. For 
a diabolic point lying exactly at the Fermi level, the same 
follows from the fact that the crossing bands are always 
on the same side of Ep in a neighbourhood of the diabolic 
point (which is always true for investigated systems, in 
which the Fermi level does not touch the borders of the 
Brillouin zone). Hence, for T — > the diabolic points 
have no singular contribution to the anomalous Hall con- 
ductivity, which we also have confirmed numerically for 
finite temperatures. It is clearly seen from © and ^ 
that the Berry curvature arises from the spin-orbit inter- 
action. This is because Hamiltonians without the spin- 
orbit coupling operator have real representations for all 
k. One can then choose u n ^ which are entirely real for 
all k and do not produce the Berry curvature. By the 
introduction of the spin-orbit coupling, the Hamiltonian 
becomes complex, causing the Berry curvature to arise. 
Yet, the diabolic points manifest themselves, when pass- 
ing the Fermi level, as kinks in the conductivity (marked 
with an arrow in Fig. [5]). 

The qualitative difference of the Berry curvature be- 
tween the models takes an effect on the anomalous Hall 
conductivity trends. The values of a xy in the k.p and 
tight-binding approaches, computed for various hole con- 
centrations p as a function of the valence band splitting 
A, are presented in Figs[5][9] The results obtained within 
the particular models remain in good agreement through- 
out the whole range of A values only for low hole con- 
centrations, p < 0.3 nm~ 3 . For higher carrier densities, 
differences in the cr xy values become significant, partic- 
ularly for small and intermediate spin splittings. Re- 
markably, we obtain a negative sign of a xy within the 
8-band and tight-binding models (Figs[6][9]) in this range: 
the higher the hole concentration, the wider the range of 
A for which the negative sign persists. This is the ef- 
fect of the Dresselhaus splitting which increases with k, 
while for increasing hole concentrations the states with 
high /c-vectors become occupied and contribute to the 
conductivity. However, strong enough spin splitting de- 
stroys the negative sign. This shows a dramatic and so 
far unnoticed influence of the Dresselhaus term on the 
AHE in hole-controlled ferromagnetic semiconductors. 

It has been suggested^ that the influence of disor- 
der on the intrinsic AHE can be phenomcnologically 
modeled by substituting one of the energy differences 
in §3§ with e rii k — e n '.k + ihT. The scattering-induced 



broadening of bands HT in (Ga,Mn)As at the localiza- 
tion boundary is presumably of the order of the Fermi 
energy Ep^lt washes out the Dresselhaus splitting and 
reduces the magnitude of its negative contribution to a xy , 
as shown in Fig. [8] However, this approach is not without 
its own problems: the energy level broadening is but a 
part of equal-rank "extrinsic" terms in the Kubo-Stfeda 
formalism^ and its magnitude is typically too large to 
treat its effect on the AHE perturbativcly. 
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FIG. 5: [color online] Anomalous Hall conductivity a xy vs 
spin splitting A for different hole concentrations p and biaxial 
strain e xx in the 6-band k.p model. The hallmark of the 
diabolic point for unstrained p = 0.5 nm -3 curve is marked 
with an arrow. Numerical errors are ca 0.5%. 
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FIG. 6: [color online] Anomalous Hall conductivity a xy vs 
spin splitting A for different hole concentrations p and biaxial 
strain e xx in the 8-band k.p model. 

The sensitivity of the AHE to the details of the band 
structure suggests that it can be influenced by the biaxial 
strain. Figures [5]|9] contain the results on the anomalous 
Hall conductivity in tensile and comprcssivcly strained 
(Ga,Mn)As samples, e xx = 1% and e xx = —1%. The a xy 
values tend to increase in the first case and decrease in 
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FIG. 7: [color online] Anomalous Hall conductivity a xy vs 
spin splitting A for different hole concentrations p and biaxial 
strain e„ in the spds* tight-binding model. 



FIG. 9: [color online] Anomalous Hall conductivity vs spin 
splitting for different hole concentrations p and biaxial strain 
e xx in the sps* tight-binding model. 
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FIG. 8: [color online] Anomalous Hall conductivity a xy vs 
spin splitting A for different hole concentrations p in the 
spds* tight-binding model. The scattering-induced energy 
level broadening KT equal to 50meV and \Ef\ is included. 



the latter in all models. As seen, small negative values 
are found already within the 6-band k.p model for the 
tensile strain. 

Additionally, we checked that despite the overall sen- 
sitivity, the effect of the temperature parameter in the 
Fermi-Dirac function on a xy for a fixed value of spin split- 
ting is negligibly small. 



ing to the experimental parameters do not fit the mea- 
sured points. The new detailed theories, which include 
the inversion asymmetry of the GaAs lattice, predict a 
negative sign of a xy for small Mn contents x. At the 
same time, all the models predict the a xy values larger 
than in the experiment for high Mn content x. They are 
nevertheless significantly lowered by the — strong in this 
regime — biaxial strain, as shown by the comparison with 
zero strain calculations for spds* tight-binding model (di- 
amonds). 
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C. Comparison to experiment 

Figure [TO] compares the theoretical and experimental 
result o 41 ' 43 on the anomalous Hall conductivity for the 
set of annealed samples with nominal Mn concentration 
x, hole concentration p and biaxial strain e xx . The calcu- 
lations of the k.p and tight-binding models corrcspond- 



FIG. 10: [color online] Reconstruction of the experimental 
anomalous Hall conductivities for nominal Mn concentration 
x, hole concentration p and strain e X x, using different theo- 
retical models. 

The much smaller than theoretical values of a xy for 
samples with high Mn content may be connected to the 
presence of Mn interstitials^ which do not form mag- 
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netic moments. The ones which survived the anneal- 
ing process, as suggested by the measured hole densities, 
compensate one substitutional Mn spin each. As a result, 
the effective Mn concentration is lower than the total Mn 
content used in our calculations, which typically leads to 
lower a xy values. This effect cannot explain the qualita- 
tive difference between the new theories and experiment 
at low Mn concentrations. However, some experimental 
data suggest that the negative a xy can be found under 
conditions indicated by the present computations^ 

Additionally, the energy levels' lifetime broadening^ 
as a part of scattering effects derived within the Kubo- 
Stfeda formalism, is taken into account in the spds* 
model (triangles). The broadening HT is taken to be the 
ratio of the total Mn concentration N = 4x/o,q and hole 
concentration p, times the magnitude of Fermi energy, 
\Ep\. As mentioned before (Fig.[8]and related text), it is 
done in a rather phenomenological way, but nevertheless 
leads to much better agreement with the experimental 
data. 



IV. SUMMARY 

We have compared four models of the (Ga,Mn)As band 
structure with regards to their impact on Curie temper- 
ature, uniaxial anisotropy and the intrinsic anomalous 
Hall effect. We considered the 8-band k.p and two tight- 
binding (spds* and sps*) parameterizations, and com- 



pared their results with the previously employed 6-band 
k.p approach. The first two quantities do not depend 
significantly on the model used, a consequence of their 
static nature. On the other hand, taking into account 
the details of the band structure beyond the six hole 
bands leads to qualitatively new results on the anomalous 
Hall effect, which is dynamic. In particular, the inversion 
asymmetry of the GaAs lattice described by the Dressel- 
haus k 3 term produces the negative anomalous conduc- 
tivity sign. Despite using the more detailed models of 
the band structure, we have not obtained the agreement 
with the experiment — indeed, moved away from it. This 
is a symptom of the intrinsic AHE theory being insuffi- 
cient to describe the observed phenomenon. Possible ad- 
ditional mechanisms which merit detailed investigation 
include scattering and localization. Their influence on 
the anomalous Hall effect in ferromagnetic semiconduc- 
tors will certainly be the subject of future studies. 
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